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TABLES  AIMD  GRAPHS  OF  THE  THEORETICAL  PEAK  PRESSURES,' 
El^RGIES,  AI®  liFuLSES  OF  SHOCK  YfAVES 
FROM  EXPLOSIYHi;  SOURCES  IN  SEA  HATER 


Abstract 


A  theory  of  the  prope.gation  of  shock  mves  from  explosive 
sooTces  ms  presented  in  NDRC  Report  A-318  (OSPD- IiBlIi),-!/  In 
that  report,  a  pair  of  ordinary  differential  equations  for  peak 
pressure  and  shock-wave  energy  as  functions  of  the  distance  from 
the  source  were  formulated  from  the  equations  of  hycFodynanics 
by  imposing  a  similarity  restraint  on  the  shape  of  the  energy¬ 
time  curve  of  the  shock  mve.  TiTO-parameter  families  of  peak 
pressure-distance  curves  are  obtained  by  the  solution  of  these 
propagation  equations.  The  parameters  are  conveniently  chosen 
as  the  initial  values  of  the  pressure  and  shock-wave  energy. 

In  the  present  report,  tables  and  graphs  of  the  two- 
parameter  families  of  curves  for  shock  ivaves  from  explosive 
sources  in  sea  mter  are  presented.  A  method  is  outlined,  for 
the  determination  of  the  parameters  from  experimental  values  of 
the  peak  pressure  and  impulse  over  a  limited  range  of  dista.nces 
from  the  source. 


I.  DESCRIPTION  OF  THE  TABLES  AND  GRAPHS  OF  THE  PROPERTIES 
OF  SHOCK  HAVES  FROM  EXPLOSIVE  SOURCES  IN  A7ATER 


1 .  Introduction 

A  theory  of  propagation  of  shock  mves  from  explosive  sources  has  been 
described  in  an  earlier  report.-^  A  pair  of  ordinary  differential  equations 
for  peak  pressure  and  shock-wave  energy  as  functions  of  distance  from  the 
source  vfere  obtained  from  the  equations  of  hydrodynamics  and  the  Hugoniot 
relations  by  imposing  a  similarity  restraint  on  the  shape  of  tlio  onergy-time 
curve  and  by  utilizing  the  second  law  of  thermodynamics  to  determine,  at  an 
arbitrary  distance,  the  distribution  of  the  initial  energy  input  from  the 
explosive  between  energy  available  for  further  propagation  and  dissipated 
energy  residual  in  the  fluid  already  traversed  by  the  shock  wavo.  Tho  theory 

—  Theory  of  the  propagation  of  shock  ivavos  from  cxplosivo  sources  in 
air  and  TJater,  by  J.  G.  Kirlovood  and  S.  R..  Brinkley,  Jr.,  I'IDRG  Report  A-318 
(0SRD-h8lii) . 
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takos  proper  account  of  the  finite  entropy  increment  in  the  fluid  resulting 
from  the  passage  of  the  shock  nave  and  permits  the  use  of  the  exact  Hugoniot 
curve  for  the  fluid. 

The  basic  equations  of  the  theory  arc 


>  (1) 


T/iiero 


§  *  -  P^=f(p), 


12Sf3 

TTT^  ^ 


12G  7  /°o"|  _ 1_ 

p-^  +  lVu/  2(1  +  ^ 


/  ^  +  2(1  - 

i(p)  =  2(1  +  g";  -  G  ^ 


G  =  1  - 


,  U  \2 
'  10  \ 


c-  =  1  £  ^ 

^  "  U  dp  ^ 

and  p  is  the  excess  peak  pressure  measured  in  units  of  a  cliaracteristic  pres¬ 
sure  Bj  at  a. .distance  Ra^  from  the  center  of  a  spherical  Ghc.rgo  of  explosive 
of  radius  a^.  B  is  the  chara.ctoristic  presstire  of  the  Tait  equation  of  state  ^ 

P=B[(y^./-l]  (2) 

and  £  is  the  exponent  of  in  the  Tait  equation..  The  qiiantitics  ^  and 

are  the  densities  of  the  Tjatcr  at  distances  R  and  infinity^  respectively;  c 
and  c^  denote  the  sound  velocity  in  vater  at  densities  £  and  respectively 
U  is  the  velocity  of  the  shock  front;  h  is  the  specific  enthalpy  increment  of 
an  element  of  fluidj  traversed  by  a  shock  -v^ve  of  peak  pressure  after  rc- 
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turn  to  pressure  p  along  its  noTr  adiabatic.  The  variable  K  is  related  to 
the  shock-TjavB  energy  £  per  gram  of  the  initial  explosive  charge  by 


t 


=  3  •^.  '1  B  „ 


(3) 


where  is  the  density  of  loading  of  the  explosive.  The  positive  impulse- 
is  given  by 


I  =  pGB., 


(U) 


where 


U  J2 
R  [G 


id 


i-^4-  +  h 


//A 


-  dR  ^ 


P  !  'O 


Tables  .of  the  functions  defined  in  Eqs.  (l)  and  (h)  are  presented  in 
an  appendix. 


2.  Description  of  the  tables  and  graphs 

The  tv/o-parancter  family  of  peak  pressure-distance  cunves  and  the  family 

of  shock-vjave  energy-distance  cui-ves  for  selected  values  of  the  constants  of 

integra/tion  were  obtained  by  numerical  integration  of  Eqs„  (l),  empl03d.ng  the 

1  / 

methods  of  Pa-rt  III  of  the  previous  report.—  Each  member  of  the  two  families 
is  characterized  by  a  value  of  p]_j  the  excess  peak  pressure  <at  the  initiai 
instant  of  time  measured  in  kilobars  and  by  a.  value  of  related  to  the 
initial  shock-^wive  energy  by 


1  Bpi  • 


(5) 


The  family  of  peak  pressure-distance  ciarvcs  is  presented  in  tabular  form  in 
Table  I.  '  In  this  and  follvvTing  tables^  a  reduced  distance  variable  equal  to 
is  employed,  T.horc  R  is  expressed  in  foot,  in  grams  per  cubic 

G  E.I  D.E  k  T  I  A  L 


contimotcr,  and  W  in  poimds.  In  order  to  facilitate  interpolation  ydth  re¬ 
spect  to  distance,  the  entries  of  the  table  are  of  the  fvinction  10  %p(^ 

2/  ® 

Yjith  p  expressed  in  pounds  per  square  inch.— 

Portions  of  the  pce.k  pressure-distance  curves  for  moderately  large  dis- 
"bcincGS  ctp  G  presented  in  a.  series  of  graphs.  Figs,  1  to  6.  Each  graph  gives 
a  family  of  peak  pressure-distance  cinrvos  for  a  fixed  value  of  the  parameter 
Qi,  plotted  as  log  Rp(n  versus  log  R(p^/ll) .  Those  graphs  arc  de- 

signed  to  aid  in  the  determination  of  the  parameters  from  an  experimental  peak 
pressure-distance  curve  in  a  manner  to  be  described  in  the  follxdng. 

The  family  of  shock-Ymive  energy-distance  curves  is  tabulated  in  Table  II, 
the  entries  of  which  are  v.dth  the  reduced  distance  defined  in  the  forego¬ 
ing  as  argument.  The  units  of  L  are  calories  pi®r  gram  of  explosive. 

The  family  _ of  positive  irapulse-distanco  curves  avas  constructed  Tdth  the 
aid  of  Eqs.  (3),  and  it  is  tabulated  in  Tabic  III.  Here,  the, function 

is  listed  iTith  the  reduced  distance  as  argioment.  The  units  of 
I  arc  pound-second  per  square  inch. 

3 .  Determination  of  the  initial  parameters  from  experimental  data 

The  graphs.  Figs,  1  to  6,  o.rc  designed  to  aid  in  the  determination  of  tho 
parameters  of  those  members  of  the  families  of  peak-pressure  curves  that  give 
the  best  fit  to  a  particular  set  of  experimental  data.  For  tlus  purpose,  the 
experimental  data  may  be  plotted  to  the  correct  scale  on  tracing  cloth  and 
compared  mth  the  theoretical  values  by  superimposing  the  tracing  on  the  pub¬ 
lished  graphs.  Unless  the  data  ai-c  of  high  accuracy,  this  procedure  mil  not 
lead  to  a  unique  determination  of  the  pe.ramctcrs  but  to  the  selection  of  sev¬ 
eral  pairs  of  values  of  the  parameters  bctwcon  which  a  choice  could  not  easily 
be  made  by  visual  methods.  The  method  of  least  squares  may  then  be  used  to 
narroiT  tho  choice  of  pa.ramotcrs.  One  computes  the  moan  of  the  squares  of  the 
vreighted  deviations  of  the  indivi.dual  experimental  values  of  the  peak  pres- 

27  = - 

—  Note  that  in  the  tables  p  is  expressed  in  pounds  per  square  inch, 
while  in  the  equations  p  m's  a  dimensionless  quantity  expressed  in  units  of  B. 


sure  and  positive  impulse  fi:om  the  theoretical  values  given  by  the  curves 
selected  by  visual  matching,  employing  the  reciprocal  of  the  root  mean  square 
deviations  of  the  experimental  points  at  a  given  distance  from  their  moan 
as  the  rnight  factor  at  that  distance.  One  then  selects  the  ;po.ramctcr  sot 

.  .  3/ 

yielding  the  smallest  value  of  this  aua.ntitj'". 

II.  TABLES  j'llD  GRAHIS  CF  THE  PEAK  PRESSURE,  ENERGY, 

AID  IliPlTLSE  OF  SHCCK  HAVES  FROH  EXPLOSIVE 
SOURCES  IN  SE/l  IViTER 


Sj^mibols  and  Units 


Symbol 

Unit 

Definition 

R 

ft 

Distance  from  the  center  of  spherical  charge  of 
explosive . 

Pe 

gm/cm^ 

Density  of  explosive. 

.W 

lb 

Height  of  charge  of  explosive. 

P 

lb /in? 

Excess  peak  pressure  of  shock  nave. 

Pi 

kilobars 

Initial  excess  peak  pressure. 

Qx 

Initial  value  of  reduced  energy  variable. 

£ 

cal/gm 

Shock-Tjave  energy. 

I 

lb-sec /in? 

Imp'alse  of  shock  vra.ve. 

Note  s 

The  first  tabulated 

corresponds  to  the 

value  of  the  distance  argument  (/)^/H)^/^=  0.136, 
surfa,cc  of  the  sphere  of  intact  explosive. 

—  See  E.  T,  I'Jhittaker “and  G,  Robinson,  The  calculus  of  observations, 
3rd  (Blackie  and  Son,  ^9hp) ,  pp.  182,  221. 
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FIG.  4.  PEAK  PRESSURE 
VERSUS  DISTANCE  CURVE 
Q.=  40. 
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APPEI^IX 


Methods  of  I-iumorical  Integra ti on 

in  the  present  section,  methods  that  are  approprie^te  for  the  nujncrical 
Integration  of  the  basic  differential  equations,  Eqs,  (l),  arc  described  and 
tables  are  given  of  the  atcciliary  functions  required  for  tlh-s  integration. 

For  shock  raves  in  rater,  it  is  convenient  to  employ  the  follovjing 
variables ; 

X  =  log  R, 

Y  =  log  P, 

Z  =  log  Q,  \  (A1 ) 

Q  =  K/Rp, 

P  =  pR. 

In  these  variables,  Eqs.  (l)  become 

d  l  _  „  ,  ^  ^  P2  ,  . 

‘  dX  " 

* 

H  =  ^(p)  +  ^  «(p), 

nhere 

^(p)  =  n(p)-1, 
m(p)  =  F(p), 
h(p)  =  F(p)~f(p). 

Tables  of  the  functions  ^(p),  m(p),  and  Yi(p)  are  presented  in  Table  IV. 

Table  V  lists  certain  auxiliary  quantities  necessary  for  the  calculation  of 
the  impulse  by  Eq.  (1|), 

A  procedure  for  the  numerical  integration  of  Eqs.  (A2)  has  been  outlined 
in  an  eo-rlier  report.—  It  is  suggested  that  the  integration  bo  started  v/ith 
the  Runge-Kutta  method,  and  an  iterative  procedure  for  continuing  the  inte¬ 
gration  is  outlined. 
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Table  IV.  Functions  for  the  integration  of  the  propagation 

equation  in  sea  water. 


p 

(Units  of  B) 

t(p) 

m(p) 

t(p) 

0 

0 

1 ,3000 

0.3000 

0.01 

0. 0011+93 

1 .1+81+3 

.1+930 

t02 

.002921. 

1.1+670 

.1+882 

.03  ' 

.001+333 

1.1+32U 

.1+839 

.oU 

.003723 

1.1+363 

.1+778 

.03 

.007073 

1 .1+207 

.1+721 

.06 

,0081+21+ 

1 .1+031+ 

.1+663 

.07 

.009771+ 

1.3871 

.1+600 

,08 

.01112 

1.3739 

.1+360 

.09 

.012U7 

1.3617 

.1+310 

0.1 

.01380 

1.31+79 

.1+1+62 

.2 

.02622 

1 .21  73 

.3981+ 

.3 

.0371+3 

1,1038 

.3376 

*h 

.01+763 

1 .0093 

.3213 

.3 

.03676 

0.9303 

.2939 

.6 

.06381 

.8333 

.2673 

.7 

,07388 

.7923 

.21+30 

.8 

.08136 

.7378 

.2268 

.9 

.08867 

.6886 

,2102 

1.0 

.09317 

.61+39 

.191+3 

.1233 

.U793 

.11+01+ 

2.0 

.11+16 

.381+1+ 

.1133 

2.3 

.1363 

.3123 

.0889 

3.0 

.1696 

.2321+ 

,0690 

3.3 

.1813 

.201+8 

.0336 

h.O 

.1919 

.1693 

.01+28 

U.3 

.2009 

.1U63 

.0366 

3.0 

.2083 

.13U8 

.031+3 

10 

.2306 

.0311+0 

,01276  - 

13 

.2678 

,02  770 

.00691+ 

20 

.281+8 

.01731 

.001+32 

23 

.311+3 

,0121+2 

.00339 

UMCLASSIFIEi: 

Table  V,  Auxiliary  table  for  blast  waves  in  sea  water. 


